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Abstract 

The study of noncommutative solitons is greatly facilitated if the field equations are integrable, 
i.e. result from a linear system. For the example of a modified but integrable U (n) sigma 
model in 2+1 dimensions we employ the dressing method to construct explicit multi-soliton 
configurations on noncommutative R 2,1 . These solutions, abelian and nonabelian, feature exact 
time-dependence for any value of the noncommutativity parameter 9 and describe various lumps 
of finite energy in relative motion. We discuss their scattering properties and prove asymptotic 
factorization for large times. 
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1 Introduction 



The past year has witnessed an explosion of activity in noncommutative field theory. The original 
motivation derived from the discovery [1] that a certain corner of string moduli space is described 
by noncommutative gauge theory. Since then, the field has been driven by the curiosity to extend 
all kinds of (quantum) field theoretical structures to the noncommutative realm [2] (for a review 
see [3]). 

Before attempting to quantize noncommutative theories, it is certainly warranted to achieve 
control over the moduli space of classical configurations. Already commutative field theories in 
1+1, 2+1, and 3+1 dimensions display a variety of nonperturbative solutions (like solitons, vortices, 
monopoles, instantons) which may be interpreted as D-branes in the string context. Turning on 
a constant magnetic NS i?-field background generalizes those branes to solitonic solutions of the 
noncommutatively deformed field theories [4, 5] (for a lecture on the subject see [6]). Since spatial 
noncommutativity requires at least two dimensions, most investigations have focused on scalar 
and gauge theories in 2+0 and 2+1 dimensions, where noncommutativity is tuned by a single 
parameter 6. One may roughly distinguish three types of results here. 

Firstly, noncommutative scalar field theories were expanded around their 6^oo limit. There, 
the potential energy dominates, and limiting static solutions are given by free projectors on sub- 
spaces in a harmonic-oscillator Fock space. Corrections due to the kinetic energy can be taken 
into account by a perturbation series in 1/6 (see [2, 7, 8, 9] and references therein). Secondly, 
explicit solutions for finite 6 were considered, mainly for abelian gauge fields interacting with scalar 
matter [2, 7 - 23]. Various properties of these configurations were investigated, including their 
stability under quantum fluctuations (see e.g. [15, 24, 25] and references therein). Thirdly, the 
adiabatic motion of solitons was studied by looking at the geometry of the moduli space of the 
static multi-soliton configurations (see e.g. [8, 9, 26] and references therein). 

In a previous paper [27] we proposed a noncommutative generalization of an integrable sigma 
model with a Wess-Zumino-Witten-type term. This theory describes the dynamics of £7(n)-valued 
scalar fields in 2+1 dimensions [28]. Its advantage lies in the fact that multi-soliton solutions to 
its field equations can be written down explicitly. This is possible because the equations of motion 
can be formulated as the compatibility conditions of some linear equations. Hence, this notion of 
integrability carries over to the noncommutative case, and powerful solution-generating techniques 
can be applied here as well. 

In the present paper we describe in detail how a new solution can be constructed from an 
old one by a noncommutative generalization of the so-called dressing approach [29, 30, 31] and 
write down a rather general class of explicit solutions, abelian as well as nonabelian, for any value 
of 6. These configurations are parametrized by arbitrary functions, and for concrete choices we 
obtain (time-dependent) multi- solitons. The latter feature finite-energy lumps in mutually relative 
motion. We discuss their functional form, large-time asymptotics, and scattering properties. 

The paper is organized as follows. We present the noncommutative extension of the modified 
sigma model, its energy functional and field equations in the following Section. Section 3 reviews 
some basics on the operator formalism of noncommutative field theory. The dressing approach is 
outlined in Section 4, where we follow our one-pole ansatz through to the general solution, which 
includes the introduction of a 'squeezing' transformation to moving-frame coordinates. Abelian and 
nonabelian static solutions including their energy, BPS bound and star-product form, are discussed 
in Section 5. We study one-soliton configurations, i.e. solutions with a single velocity parameter, in 
Section 6, derive their energies and give examples and limits. Section 7 finally exemplifies the multi- 
soliton construction on the two-soliton case (with relative motion), abelian as well as nonabelian. 
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The asymptotic behavior at high speed or large times is derived, which aids the evaluation of the 
energy. In three appendices, we collect formulae about coherent and squeezed states and present 
two different proofs of the large-time factorization of multi-soliton configurations. 

2 Modified sigma model in 2+1 dimensions 

Definitions and notation. As has been known for some time, nonlinear sigma models in 2+1 
dimensions may be Lorentz-invariant or integrable but not both. In this paper we choose the second 
property and investigate the noncommutative extension of a modified U (n) sigma model (so as to 
be integrable) introduced by Ward [28]. 

Classical field theory on noncommutative spaces may be realized in a star-product formulation 
or in an operator formalism. The first approach is closer to the commutative field theory: It is 
obtained by simply deforming the ordinary product of classical fields (or their components) to the 
noncommutative star product 

(f*9)(x) = f(x)exp{^d a 9 ab d b }g(x) , (2.1) 

with a constant antisymmetric tensor 9 ab , where a, b, . . . = 0, 1, 2. Specializing to R 1 ' 2 , we shall use 
(real) coordinates (x a ) = (t,x,y) in which the Minkowskian metric reads (r) a b) = diag(— 1, +1, +1). 
For later use we introduce two convenient coordinate combinations, namely 

u := \{t + y) , v := \{t - y) , d u = d t + d y , d v = d t - d y (2.2) 

and 

z := x + iy , z := x - iy , d z = \{d x -id y ) , d- = \{d x + \d y ) . (2.3) 

Since the time coordinate remains commutative, the only non- vanishing component of the noncom- 
mutativity tensor is 

gxy = _gyx = . Q > Q ^ Q zz = _ Q zz = _ 2[0 ( 2 ^ 



Action and energy. The noncommutative U (n) sigma model may be obtained by a reduction of 
the Nair-Schiff sigma-model-type action [32, 33] from four to three dimensions, 

S = ~ Jdtdxdyri^tr^da®- 1 * d 6 $) 

- ^JdtdxdyJ^ dpv x e x ^ a tr ★ * ★ d u $> ★ ★ d a ^j , (2.5) 

where Greek indices include the extra coordinate p, and e XtJ,ua denotes the totally antisymmet- 
ric tensor in IR 4 . The field <fr(t,x,y) is group-valued, & = < I>~ 1 , with an extension <&(t,x,y,p) 
interpolating between 

®(t,x, y, 0) = const and x,y, 1) = $(t,x,y) , (2.6) 



and 'tr' implies the trace over the U(n) group space. Finally, (v\) = (v c , 0) is a constant vector in 
(extended) space-time. For (v c ) = (0, 0, 0) one obtains the standard (Lorentz-invariant) model. Fol- 
lowing Ward [28], we choose (v c ) = (0, 1,0) spacelike, which yields a modified but integrable sigma 
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model. Although this explicitly breaks the Lorentz group of SO(l, 2) to the GL(1, R) generated by 
the boost in y direction, it leaves unmodified the conserved energy functional 

E = ^Jdxdy tr (d t &*d t & + d x &*d x § + d y ^*d y ^j . (2.7) 

Field equations. The noncommutative sigma-model equation of motion following from (2.5) reads 

( V ab + v c e cab )d a ($- 1 *d b <l>) = , (2.8) 

where e abc is the alternating tensor with e 012 =l. In our coordinates (2.2) it is written more concisely 
as 

a*^- 1 *^)-^^ -1 *^) = . (2.9) 

This Yang-type equation [34] can be transformed into a Leznov-type equation [35], 

d x (j) — d u d v <f> + d v (f> * d x (j) — d x (j> * d v (j> = , (2-10) 
for an algebra-valued field (f>(t, x, y) G u(n) defined via 

d x( j) ■= fc" 1 *^ =: A and d v <f) := * _1 *S X $ =: B . (2.11) 

The equation (2.9) actually arises from the Bogomolnyi equations for the 2+1 dimensional 
Yang-Mills-Higgs system (A^ , ip) , 

habcF bc = d a <p + A a *^- V *A a . (2.12) 

Indeed, after choosing the gauge A v = = A x +(p and solving one equation by putting 

A u = $ -1 *0„$ and A x - <p = $ _1 *a a: $ , (2.13) 

the Bogomolnyi equations (2.12) get reduced to (2.9). 

3 Operator formalism 

Fock space. The nonlocality of the star product renders explicit computations cumbersome. 
We therefore pass to the operator formalism, which trades the star product for operator-valued 
coordinates x M satisfying [f*,^] = \9^ v '. The noncommutative coordinates for R 1,2 are (t,x,y) 
subject to 

[t,x] = [t,y] = , [x,y] = W => [z,% = 29 . (3.1) 

The latter equation suggests the introduction of (properly normalized) creation and annihilation 
operators, 

z and = —= z so that [a, at] = 1 . (3.2) 



'IB V26 

They act on a harmonic-oscillator Fock space TL with an orthonormal basis {\n), n = 0, 1,2, . . .} 
such that 

a)a\n) =: N \n) = n\n) , a\n) = y/n\n— 1) , a'\n) = y/n+1 \n+l) . (3.3) 
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Moyal-Weyl map. Any function f(t,z,z) can be related to an operator-valued function f(t) = 
F(t,a,(v) acting in TL, with the help of the Moyal-Weyl map (see e.g. [36, 24]) 

f{t, z,z) — F(t,a,«t) = -J^dzdzf^z^e-^"-* 

= Weyl-ordered f(t, V29a, VW<J) . (3.4) 
The inverse transformation recovers the c-number function, 

F(t,a,a<) — . f(t,z,z) = 2iT0 f^S ^{f((, a, a') e i|p(V2»«-.)+p(V2»«t_ z)] j 

where 'Tr' signifies the trace over the Fock space H, and F* is obtained from F by replacing ordinary 
with star products. Under the Moyal-Weyl map, we have 

/*<? — /$ and [dxdyf = 2vr#Tr/ = 2vr^(n|/|n) . (3.6) 

n>0 

The operator formulation turns spatial derivatives into commutators, 

d x f — ^[y,/] =: 0*/ and 0„/ — ~[x,f] =: 0„/ , (3.7) 

so that 

0*/ dj = [at,/] and ^/ £_/=_!=[«,/] . (3.8) 



(3.9) 



The basic examples for the relation between / and / are 
/ : 1 2 z zz— 6* 

/: 1 V26a ^29 a) 29N 26(N+±) 26(N+1) 

For more complicated functions it helps to remember that 

z-kz = zl + 9 and z*z = zz — 9 , (3.10) 

because the composition law (3.6) allows one to employ the star product (on the functional side) 
instead of the Weyl ordering (on the operator side). For notational simplicity we will from now on 
omit the hats over the operators except when confusion may arise. 

4 Dressing approach and explicit solutions 

The payoff for considering an integrable model is the availability of powerful techniques for con- 
structing solutions to the equation of motion. One of these tools is the so-called 'dressing method', 
which was invented to generate solutions for commutative integrable systems [29, 30, 31] and is 
easily extended to the noncommutative setup [27]. Let us briefly present this method (already in 
the noncommutative context) before applying it to the modified sigma model. 

Linear system. We consider the two linear equations 

{(d x -d u )ip = Ail) and (0 V - d x )1> = , (4.1) 
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which can be obtained from the Lax pair for the (noncommutative) self-dual Yang-Mills equations 
in M 2 ' 2 [37, 38] by gauge-fixing and imposing the condition d^ip = 0. Here, ip depends on (t, x, y, Q 
or, equivalently, on (x,u,v,Q and is an nxn matrix whose elements act as operators in the Fock 
space TC. The matrices A and B are of the same type as ip but do not depend on £. The spectral 
parameter £ lies in the extended complex plane. The matrix ip is subject to the following reality 
condition [28]: 

il>{t,x,y,Q [V(i,x,y,C)]t = 1 , (4.2) 
where 'f is hermitian conjugation. We also impose on tp the standard asymptotic conditions [38] 

^(t,x,y,C^oo) = 1 + fVM,!/) + 0(C 2 ) , (4.3) 
4>(t,x,y,(^0) = $-\t,x,y) + O(0 . (4.4) 

The compatibility conditions for the linear system of differential equations (4.1) read 

d x B-d v A = , (4.5) 

d x A - 8 U B - [A, B] = . (4.6) 

We have already encountered 'solutions' of these equations: Expressing A and B in terms of <f> like 
in (2.11), 

A = d x <p = $ -1 d u $ and B = d v (f> = d x $ , (4.7) 

solves the first equation and turns the second into (the operator form of) our Leznov-type equa- 
tion (2.10), 

d 2 x 4>-d u d v 4>- [d^, d v 4>] = o . (4.8) 

Alternatively, the ansatz employing $ in (4.7) fulfils the second equation and transforms the first 
one into (the operator version of) our Yang-type equation (2.9), 

d x ($~ 1 d x $) - d v ($~ 1 a M $) = . (4.9) 

Inserting these two parametrizations of A and B into the linear system (4.1) we immediately verify 
(4.3) and (4.4), confirming that the identical notation for A, B, 4>, and $ in Section 2 and in the 
present one was justified. 

Dressing approach and ansatz. Having identified auxiliary linear first-order differential equa- 
tions pertaining to our second-order nonlinear equation, we set out to solve the former. Note that 
the knowledge of ip yields (j) an d by way of (4.3) and (4.4), respectively, and thus A and B 
via (4.7) or directly from 

-W,x,y,Q(Cd x -d u )[il>(t,x,y,d]t = A(t,x,y) , (4.10) 
-i;(t,x,y,C)(Cd v -d x )m,x,y,C^ = B(t,x,y) . (4.11) 

The dressing method is a recursive procedure generating a new solution from an old one. Let us 
have a solution ipo(t, x,y, Q of the linear equations (4.1) for a given solution (Aq,Bq) of the field 
equations (4.5) and (4.6). Then we can look for a new solution tp in the form 

s m 

4>(t,x,y,() = x(t,x,y,Qifo(t,x,y,Q with x = 1 + . kot )a , (4.12) 

a=l fc=l ^ ^ k > 

where the Hk(t,x,y) are complex functions and the nxn matrices Rk a (t, x, y) are independent of (. 
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In this paper we shall not discuss this approach in full generality (for more details see e.g. [29, 
30, 31]). We consider the vacuum seed solution Aq = Bq = 0, ipo = 1 and a restricted ansatz 
containing only first-order poles in (: 

i> = X^Po = 1 + , (4-13) 

P =i <> 

where the \i v are complex constants with Im/i^ < 0. Moreover, following [31], we take the matri- 
ces R k to be of the form 

m 

R k = Yj T e r * T l ( 4 - 14 ) 

t=\ 

where the T k (t,x,y) are nxr matrices and F ik (t,x,y) are rxr matrices with some r>l. In the 
abelian case (n=l) one may think of T k as a time-dependent row vector (\z\,t), \zl,t), . . . , \z r k ,t)) 
of (ket) states in H. 

We must make sure to satisfy the reality condition (4.2) as well as our linear equations in the 
form (4.10) and (4.11). In particular, the poles at £ = fi k on the l.h.s. of these equations have 
to be removable since the r.h.s. are independent of £. Inserting the ansatz (4.13) with (4.14) and 
putting to zero the corresponding residues, we learn from (4.2) that 

m „ 

(l-£-^K = (4.15) 
while from (4.10) and (4.11) we obtain the differential equations 

{Aotb) = -v>(C)«(C)t _ (i-Err-) 1 ^! = ( 4 - 16 ) 

Q-^Vk ^ -y I^P l^k ' 

where 

L k := Jx k d x - d u or L k := fl k d v - d x . (4.17) 
The algebraic conditions (4.15) imply that the T £p invert the matrices 

-. m 

r pk = ^TlT k , i.e. Vr% = ^ . (4.18) 

Hp-Vk P p[ 

Finally, by substituting (4.13) and (4.14) into the formulae (4.4) and (4.3) we solve the equations 
of motion (4.9) and (4.8) by 

m -■ m 

$-! = $t = i _ ^ —T e T ek Tl and ^ = Yj Te^tf (4.19) 

k,t=\ ^ k k,i=\ 

for every solution to (4.15) and (4.16). 

Moving-frame coordinates. At this point it is a good idea to introduce the co-moving coordi- 
nates 

w k := u k [x + p, k u + fL^v] = v k [x+ \(jl k - Ji k V )y + \{jlk + V k 1 )t] 

w k := v k [x + fi k u + ^v] = v k [x+ \{n k - fj, k r )y + |(// fe + /i^ 1 )t] , (4.20) 



6 



where the v k are functions of \x k via 

4i /',, - //, ' - 2i 1 1 2 



(4.21) 



In terms of the moving-frame coordinates the linear operators defined in (4.17) become 

Lk = v k (jl k - ^k)0 Wk or L k = VkV^ifih- Vk)dw k , (4.22) 
respectively, so that we have only one equation (4.16) for each pole. Since 

[w k ,w k ] = \9 v k v k (fi k - p, k - /i^ 1 + fl^ 1 ) = 29 (4.23) 
we are led to the definition of co-moving creation and annihilation operators 



c fc = ~^w k and c\ = —=^w k so that [c k ,c\] = 1 . (4.24) 



1 ~ A t 1 

= Wk and c! = —= 
29 k V26 

The static case (w k =z) is recovered for [i k = — i. Coordinate derivatives are represented in the 
standard fashion as [36, 24] 

V26d Wk — -[c\,.] , V26d Wk — [ Cfc ,.] . (4.25) 

The number operator N k := c k c k is diagonal in a co-moving oscillator basis {\n) k ,n=0, 1,2, . . .}. 

It is essential to relate the co-moving oscillators to the static one. Expressing w k by z via (4.20) 
and (4.21) and using (3.2) one gets 

c fc = (coshrfc)a- (e^sinhTfc)^ - (3 k t = U k (t)aUl(t) , (4.26) 
which is an inhomogeneous SU(1, 1) transformation mediated by the operator [39] 

U k (t) = e ^ a k a12 ~^ka 2 e (/3 k a*> - j3 k a)t _ (4.27) 

Here, we have introduced the complex quantities 

a k = e^Tfc and (3 k = -\v k (j2 k +j2 k l )/V29 , (4.28) 

with the relations 

£fc := e 1,?fc tanhTfc = ^ fc ~^fc " 2l and ^ = cos h Tfc - e 1,?fe sinh r fc , (4.29) 

Hk~H k + 21 

so that [i k may be expressed in terms of £ k via fi k = i^y^p-. Clearly, we have |n)fc = U k (t)\n), and 
all co-moving oscillators are unitarily equivalent to the static one. The unitary transformation (4.26) 
is known as 'squeezing' (by a k ) plus 'shifting' (by [3 k t). 

Solutions. After this diversion, the remaining equations of motion (4.16) take the form 

m p 

(i-ETTzV)^ = > ( 43 °) 

v p=1 H'k / 
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which means that c k T k lies in the kernel of the parenthetical expression. Obviously, a sufficient 
condition for a solution is 

c k T k = T k Z k (4.31) 
with some rxr matrix Z k . We briefly elaborate on two important special cases. 
Firstly, if n>2 and 

Z k = lc fc then [c k ,T k ] = , (4.32) 

which can be interpreted as a holomorphicity condition on T k . Hence, any nxr matrix T k whose 
entries are arbitrary functions of c k (but independent of cl) provides a solution ib of the linear 
system (4.1), after inserting it into (4.18), (4.14), and (4.13). The complex moduli \i k turn out to 
parametrize the velocities v k of individual asymptotic lumps of energy. The shape and location of 
these lumps is encoded in further moduli hiding in the functions T k (c k ). 

Secondly, the abelian case (n=l) deserves some attention. Here, the entries of the row vector 
T k = {\z\,t), \z\, t),..., \z k ,t)) are nonfunctions of c k . If we take 

Z k = diag(4,^,...,4) with 4eC then c k \z{, t) = z\ \z\, *) , (4.33) 

which explains our labeling of the kets and qualifies them as coherent states based on the co-moving 
ground state |0)fc, 

|4, t) = e ^ c I-^ Cfc |0) fc = e 44-4 c * U k (t) |0) = U k (t)e z ^~^ a \0) . (4.34) 

In total, the abelian solutions \z k ,t) are nothing but squeezed states. The squeezing parameter 
is a k while the (complex) shift parameter is z\ + (3 k t. In the star-product picture, we find m 
groups of r asymptotic lumps of energy. The squeezing parametrizes the (common) deviation of 
the lumps of type k from spherical shape, while the shifts yield the position of each lump in the 
noncommutative plane. Clearly, all lumps in one group move with equal constant velocity but are 
arbitrarily separated in the plane. This discussion extends naturally to the nonabelian case. 

We remark that our solutions have a four-dimensional perspective. Recall that the linear sys- 
tem (4.1) can be produced from the linear system of the self-dual Yang-Mills equations on noncom- 
mutative M 2 ' 2 . Therefore, the dressing approach may be employed in this more general situation 
to create BPS-type solutions of the functional form (4.19) to the self-dual Yang- Mills equations 
in 2+2 dimensions. Hence, our soliton solutions of (4.9) in 2+1 dimensions can be obtained from 
these infinite-action BPS-type solutions by dimensional reduction. In contrast, BPS-type solutions 
of the form (4.19) do not exist on noncommutative M 4 ' because the reality condition (4.2) for the 
matrix ib then involves £ — ► — ~ instead of £ — ► £. However, in four Euclidean dimensions a modified 
ADHM approach allows one to construct noncommutative instanton solutions the first examples of 
which were given by Nekrasov and Schwarz [40]. 

5 Static solutions 

Projectors. For comparison with earlier work, let us consider the static case, m = 1 and fi = — i. 
Staying with our restricted ansatz of first-order poles only, the expression (4.13) then simplifies to 

ib = 1 + =: 1 P so that <J> t = l- 2P = P J _-P . (5.1) 

C+i C+i 

The reality condition (4.2) directly yields 

P ] = P and P 2 = P , (5.2) 
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qualifying P as a hermitian projector. Indeed, (4.2) degenerates to (l—P)P = 0, while (4.16) 
simplifies to 

(l-P)LP = =► F := (l-P)aP = . (5.3) 

We emphasize that this equation also characterizes the BPS subsector of general noncommutative 
scalar field theories [8, 9]. 

It is illuminating to specialize also our ansatz (4.14) for the residue, with (4.18), to the static 
situation: 

T = ^ P = — = T — Tt = T-Lpt . (5.4) 

TtT -2i -2i TtT v ; 

Hence, the projector P is parametrized by an nxr matrix T, where, in the nonabelian case, r<n 
can be identified with the rank of the projector in the U(n) group space (but not in 7i). The 
equation (4.30) then simplifies to 

(l-P)oT = , (5.5) 
which means that oT must lie in the kernel of 1—P. 

Energy and BPS bound. Clearly, within our ansatz static configurations are real, <3?t = <!>. The 
energy (2.7) of solutions to (5.3) reduces to 

E = 4:irOTr(d z &&z®) 

= 87rTr ([a<,P] [P, a]) = 87rTr (at Pa - Pa 1 " a + 2a 1 " f) 

= SvrTrfVPa-Pafa) = 8vrTrP , (5.6) 

where the last equality holds only for projectors of finite rank in H. 

The topological charge Q derives from a BPS argument. The energy of a general static config- 
uration (within our ansatz but not necessarily a solution of F=0) obeys 

±E = Tr(y,P][P,a]) = Tr (p [a f , P] [a , P] - P [a , P] [a f , P]) + Tr (Vf + P F f ) 
> Tr (P [at, P] [a, P]-P [a, P] [at, P]) 

= §0Tr ($^$#2/$ - $S y *S x *) == Q ■ (5-7) 
Alternatively, using G:= (1-P) at Pin 

E = Tr ([at, P][P, a]) = -Tr (p [a*, P] [a , P] - P [a , P] [at, P]) + Tr (&G + GG^) (5.8) 
yields, along the same line, E > — 8ir Q. In combination, we get 

E > 8vr|Q| , with E = 8tt\Q\ iff P = or G = . (5.9) 

The P=0 configurations carry positive topological charge (solitons) whereas the G=0 solutions 
come with negative values of Q (antisolitons). 

It is instructive to imbed the BPS sector into the complete sigma-model configuration space. 
For static configurations, the sigma-model equation (4.9) becomes 

d g (&dt$) + d,(&d g *) = [[a,[at,P]],P] = . (5.10) 
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We emphasize that this is the equation of motion for the standard noncommutative Euclidean 
two-dimensional grassmannian sigma model, as was clarified previously [27]. A short calculation 
shows that solitons (F=0) as well as antisolitons (G=0) satisfy this equation. To summarize, our 
choice of ansatz for ip leads to the construction of solitons. For describing antisolitons it suffices to 
switch from n\~ to Jik- 

Abelian static solutions. Abelian solutions (n=l) are now generated at ease. Since fi = — i 
implies a=0 and (3=0, the squeezing disappears and r lumps are simply sitting at fixed positions 
z l , with i=l, . . . ,r. Hence, we have 

T = (|z 1 ),|z 2 ),...,|* r )) with \z l ) = e* <aWa |0) . (5.11) 

This is obvious because w — > z in the static case. The coherent states \z % ) are normalized to unity 
but are in general not orthogonal. One reads off the rank r projector 

r r—1 

P = Y.\ zi )^ z '\ z ')Yi-(* j \ = wX>><n|Wt , (5.12) 

i,j=l n=0 

where the last equality involving a unitary operator U was shown in [9]. Clearly, the row vector 
T = (|0), |1), . . . , |r— 1)) directly solves (5.5), and P projects onto the space spanned by the r lowest 
oscillator states. Obviously, these solutions have E = 8irr. 

Nonabelian static solutions. Another special case are nonabelian solutions with r=l. They are 
provided by column vectors T whose entries are functions of a only. If we choose polynomials of 
at maximal degree q, the energy of the configuration will be E = 8nq, independently of the gauge 
group. To illustrate this fact we sketch the following U{2) example: 



T = 



_ ^ -= * \ 

XX+Nl/(N-q)\ XX+N\/(N-q)\ 

,<? _= A a i i a \ q i 

' \X+Nl/(N-q)\ XX+N\/(N-q)\ j 



(5.13) 



which is of infinite rank in TC. Still, the energy (5.6) is readily computed with the result of 8nq. 

Inverse Moyal-Weyl map. Having constructed noncommutative solitons in operator language, 
it is natural to wonder how these configurations look like in the star-product formulation. The 
back translation to functions &*(t,x,y) solving the noncommutative field equations (2.9) is easily 
accomplished by the inverse Moyal-Weyl map (3.5). We take a brief look at r=l solutions. 

The abelian case has been considered many times in the literature [2, 15, 16, 6]. Here one has 
$ = 1 - 2|0)(0| — ► = 1 - 4e~ zl/e , (5.14) 

which lacks a nontrivial #^0 or 9^oo limit. 

For a nonabelian example we consider the static U(2) solution <E> = 1 — 2P with P given 
by (5.13), for degree q=l. With the methods outlined in Section 3, we find 

xx+n xx+N a \ ^ I fy+zz-e ^ j^+^wp z 

I — > ®* = ( 5 - 15 ) 

n x i t / \ 2z -ri+ zz - 26 „ ~ I ii-zz-e 

AA+7V AA+]V / V ZZ (^+zz~e)^ + 7 -f+zz+e 



with 7 = V29X. One may check that it indeed fulfils (2.9). 
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It is instructive to display the 9 — > and 9 — > oo limits of <!>*, while keeping 7 fixed: 



77—22 —272 
77+22 77+22 



+1 

$*(0-O) - I _ "_ I and $*(0^oo) ^ I I . (5.16) 



-227 , 77—22 

77+22 ' 77+22/ 



-1. 



Other than the abelian solution, $+(z, z, 9) is a rational function with a nontrivial commutative 
limit. 

6 One-soliton configurations 

For a time-dependent configuration with m=l, our first-order pole ansatz (4.13) simplifies to 
ij> = 1 + — — =: 1 + ^^P with P = T-J-Tt , (6.1) 
where P and T satisfy 

(l-P)cP = and cT = T Z . (6.2) 

Here, 

c = (cosh r) a- (e k sinh r) a f - (3t = U{t)aU ] {t) (6.3) 

is the creation operator in the moving frame. 

Comparing the formulae of Section 4 with (6.1)-(6.3) we see that the operators c, c' and there- 
fore the matrix T and the projector P in (6.2) can be expressed in terms of the corresponding 
static objects (see Section 5). This is accomplished by means of the inhomogeneous SU(1, 1) 
transformation (6.3), which on the coordinates reads 

(w\ I coshr — e^sinh-A I z\ 1 fv (fi + /i _1 )\ 
w) V-e-^sinhr coshr ) \z) + 2 [v ( M + ^ ) 1 

and is in fact a symplectic coordinate transformation, preserving dz A dz. Even though the inho- 
mogeneous part of (6.4) is linear in time, this does not mean that the solution 

$ = $- = 1-pP with p=l-p,//j, (6.5) 

is simply obtained by subjecting the corresponding static solution to such a 'symplectic boost'. 
To see this, consider the moving frame with the coordinates w,w and t'=t and the related change 
of derivatives, 

d x = vdw + vdyj , 

d t = ±v {fi + fi' 1 ) d w + \v {pi + iT 1 ) d w + d t , , 

d y = \v{ji - JT 1 ) d w + \v (p - fi" 1 ) duj . (6.6) 

In the moving frame our solution (6.5) will be static, i.e. dt'<& = 0. However, the field equation (4.9) 
is not invariant w.r.t. the coordinate transformation (6.4). Therefore, the solution (6.5), which is 
static in the moving frame, has a functional form different from that of the static solution (5.1) - 
the coefficient of the projector is altered from 2 to p. Due to this non-invariance, the coordinate 
frame (z,z,t) is distinguished because in it the field equation takes the simplest form. 
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The constant velocity v in the xy plane is easily derived by writing down the map (6.4) in real 
coordinates (w=:x'+iy'), 1 



x' 



(cosh r — cos i? sinh r — sin t? sinh r ' 
— sin $ sinh r cosh r + cos d sinh r/ L\y, 



(6.7) 



where we find 

"(S • ' (6 - 8) 

The velocity dependence of the energy 

E = neTr(g ab d a $td b $i>) with Qf 6 ) = diag (+1, +1, +1) (6.9) 
is now obtained straightforwardly by employing derivatives w.r.t. w and W via (6.6), which yields 

Inserting this bilinear form and (6.5) into the energy functional (6.9), we are left with [27] 

= Trei6f(v)Tr(d w Pd w P^ 
= 87r/(i/)Tr([ct,P][P,c]) 

= 87Tf(v)q = f(v)E[^} , (6.11) 

because the final trace is, by unitary equivalence, independent of the motion. We have introduced 
the 'velocity factor' 



_i 1 Vl — v 2 

f(v) := jf- pp vv (y/JIfi + y/JIfi ) 2 = — = cosh rj sin ip = — 5- . (6.12) 

vv 1 — v~ 

It is noteworthy that for motion in the y direction (<p=^) the 'velocity factor' cosh 77 agrees with 
the relativistic contraction factor l/y/l— v 2 > 1, but it is smaller than one for motion with \v x \ > 
\v y \y/l— v'y, and the energy may be made arbitrarily close to zero by a large 'boost', since v 2 — > 1 
is equivalent toc^^0or?7^±oo. 

Examples. We can immediately write down moving versions of our solutions (5.14) and (5.15) 
from the previous section. Since |0)' := U(t)\0) denotes the moving ground state (i.e. squeezed and 
shifted), the r=l abelian soliton reads 

$ = 1 - pU(t)\0)(0\U(tf => = 1 - 2pe~ wW/e 

= 1 - 2 pe-F-^ TATA F-* t V6 j ( 6 .i3) 



x For geometric visualization, the particular 5 I L(2,R) transformation matrix appearing here may be factorized to 
cos 8 sni(5\ /coshr sinh-A /cos 8 — sin<5N 



sin 8 cos 51 \ sinh r cosh r ) \ sin 5 cos 8 



where 25 = -f - 1?. 
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where A is the 'squeezing matrix' appearing in (6.7). In a similar vein, the r=l nonabelian solu- 
tion (5.15) gets generalized to 

/ _aa _A t \ / _aa _a t \ 

. , / AA+ata AA+ata \ „ , s+ / AA+ctc AA+ctc \ 

$ = 1 - pU(t) _ U{tf = 1 - p _ , (6.14) 

\ AA+ata AA+ata / \ AA+ctc AA+ctc / 

which leads to 

( 77_ 7 2 +w ^~ 26> , TlA 
77+m>m>— # (77+uiui— 0) 2 \ 

, (6.15) 
77+ ww- 26 — _«iin+fl / 

(77+touJ— 6) 2 ' ff+ww+9 / 

with the time dependence hiding in 

w = w(t,x,y) = (cosh r — e 11? sinhr) (x — v x t) + i (cosh r + e w sinhr) (y — v y t) , (6.16) 
as expected. We note that in the distinguished frame the large-time limit (for v^O) is 

lim = 1 - pU with n = [ ? ? ) . (6.17) 



Although a one-soliton configuration with topological charge q>l may feature several separated 
lumps in its energy density it does not deserve to be termed a 'multi-soliton', because those lumps 
do not display relative motion. True multi-solitons reduce to a collection of solitons only in the 
asymptotic regime (\t\ — > 00) or for large relative speed (v Te \ — ► 00); they will be investigated next. 

7 Multi-soliton configurations 

Two-soliton configurations. In Section 4 we gave an explicit recipe how to construct multi- 
soliton configurations 

TO 

d>t = 1 - J2 R */Vk (7.1) 

fe=l 

via (4.14), (4.18), and (4.30). It is known from the nonabelian commutative situation [28] that 
for polynomial functions T^(cfc) such solutions have finite energy and describe m different lumps 
of energy moving in the xy plane with distinct velocities vj. parametrized by p k =e rik ~ ltfk , in other 
words: multi-solitons. Since the explicit form of such configurations is more complicated than 
that of one-soliton solutions, we illustrate the generic features on the simplest examples, namely 
two-soliton configurations with r=l. Thus, in the following, m = 2 and k = 1, 2. 

We begin by solving T £p in terms of via (4.18), taking into account the noncommutativity: 

+P+1.1 — r^r^i^i] 1 , r 12 = — [fn — r^r^i^i] 1 ri2r 2 2 1 > 

-[f2 2 -r 2 ifn 1 ri2]" 1 f2ifn 1 > T 22 = +F22 - f 21 f n}' 1 , (7.2) 
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which in terms of T\ and T 2 reads 



pll 


= + 


" 1 


t\t x - 


M22 










pl2 




' 1 


t\t x - 


M22 






Mil 


M12M21 


p21 




■ 1 


t\t 2 - 


Mil 






M22 




M21M12 


p22 


= + 


" 1 


t\t 2 - 


Mil 






-M22 


M21M12 



1 



T±T 2 —l j- 2 

rfTa-p-rtri 



T{Ti 
1 



M22 T t T2 1 

M12 T 2 T T 2 
Mil „t^r 1 

J 2 J 1 + 

M21 T/Ti 



t\t 2 



(7.3) 



with (i k e 



Mfc — W- For the full solution, we then arrive at 
= 1 - Tir 11 ^/^! - Tir 12 T 2 V//2 - T 2 T 21 Tl/^ - T 2 T 22 T^/fi 2 

= 1 _ f^T 1 [Tl(l-aP 2 )T 1 ]- 1 Tl + ^i (J T 1 [T 1 t (l-aP 2 )T 1 ]- 1 T 1 t P 2 
Mi M2 

+ f^aT 2 [T^(l-aP 1 )T 2 ]- 1 T^P 1 - ^ T 2 [T^l-aPi)^]" 1 T 2 f 
Mi " " M2 

1 1 [Mil p M21 DD \ 1 fM22 p M12 DD \ 
= 1 - - 5-5- 4 Pi cr PiP 2 f - — — \ P 2 a P 2 Pi \ 



(7.4) 
(7.5) 



where a := 



G R, the T fc must be chosen to satisfy (4.30), and P k = Tfc-^-T* as usual. 



The abelian rank-one two-soliton is included here: One simply writes 



\zk,t) 



Zkcl-ZkCk 



Zkcl-ZkCk 



U k (t)\0) = U k (t)e z ^-^ a \0) 



(7.6) 



Thus, noncommutative abelian fields look like commutative U(oo) fields. Our kets are normalized 
to unity, and therefore P k = \zk,t)(zk,t\. In the bra-ket formalism our solution takes the form 



1 1 / Mil 1 ,\ / , 1 M21 1 , w , 

1 - — tk \ \zi,t){zi,t\ - as\zi,t){z 2 ,t 

l-a \s\ z I mi fi 2 



- —as\z 2 ,t)(z 1 ,t\ + —\z 2 ,t)(z 2 ,t\\ , (7.7) 
Mi M2 J 

where the overlap (zi,t\z 2 ,t) was denoted by s. 

The simplest nonabelian example is a U (2) configuration with T\ = ( \ ) and T 2 = (\) = UT\U\ 
i.e. taking Mi = — i and [i 2 =: M- We refrain from writing down the lengthy explicit expression 
which results from inserting this into (7.3) and (7.4) or from plugging the projectors P k = T k ^ Tj, 

into (7.5). 

High-speed asymptotics. The two-soliton configuration (7.4) does not factorize into two one- 
soliton solutions, except in limiting situations, such as large relative speed or large time. Let us 
first consider the limit where the second lump is moving with the velocity of light, 



<P 



or 7T 



Im [i 2 







\v 2 \ 



1 



(7.8) 
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while the velocity of the first lump is subluminal. In this limit, we have a — > and therefore 

*t _> i _ + ^ a P x P 2 + ^aP 2 P, - ^P 2 . (7.9) 

Mi M2 Mi M2 

Define also p^ ■= Pkk/pk = ^—pk/pk- Clearly, Im p 2 — ► implies that 

$t _> i _ piPl , (7.10) 
except for v 2 — > (0, ±1) meaning /i 2 — ► or oo, where 

$ f (l-piPi)(l-/9 2 P2) or $t ^ (l-p 2 P 2 )(l- Pl P 1 ) , (7.11) 

respectively. 

Large-time asymptotics. Now we investigate the behavior of (7.4) for large (positive and neg- 
ative) time. For convenience, we take the first lump to be at rest while the second one moves, 
i.e. 

Mi = -i and M2 = : M =► Pi = 2 , P2 =■ P , a = 1 - ^ — • ( 7 - 12 ) 

(M-i)(M + i) 

For simplicity, we consider r=l, i.e. P& is a column vector and P& of (matrix) rank one. After some 
tedious algebra (see Appendix B) one learns that <&' factorizes in the large-time limit: 

$t ^ (1-2 Pi) (l-pn 2 ) , (7.13) 

where II 2 := limi t i_ >0O P 2 is a coordinate-independent hermitian projector, and 

Pi = [an 2 + (i-n 2 )]Ti n : — rf [an 2 + (i-n 2 )] with a = ^— - (7.u) 

T{ (l-o-II 2 ) T\ P~ l 
is a new projector describing a one-soliton configuration. 

We have explicitly checked this behavior for the simple U (2) example mentioned above and 
found indeed 

$t(z,y,t->±oo) = (1-2 Pi) (l- P U 2 ) + Oir 1 ) with Pi = fi^Lf/ , (7.15) 

Pi Pi 

where 

' 1 A , A) 0' 



a ay and H 2 = ^ ^ . (7.16) 

Let us now consider the situation within the frame moving with the second lump. In this 
frame the second lump is static. By the same arguments as above, the solution <J>t behaves in the 
large-time limit as 

$t ^ (1-2IL) (l-pP 2 ) , (7.17) 

where IIi is a coordinate-independent hermitian projector, and the projector P 2 describes the 
second lump in the moving frame. It is obvious how these statements generalize to m-soliton 
configurations. If we sit in the frame moving with the Ith lump, then the large-time limit of 3>t is 
a product of (1 — pePg) and constant unitary matrices, 

$t ^ ^ _ piHl ) (1 _ p2 u 2 ) . . . (1 - W _iILj_i) (1 - Pt P t ) (1 - pt+iSh+i) ■ ■ ■ (1 - Pm^m) ■ (7.18) 
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The proof is outlined in Appendix C. 

A lesson from (7.18) concerns the scattering properties of our multi-soliton solutions. Since the 
limits t — > — oo and t — > +00 in (7.18) are identical, every lump escapes completely unharmed from 
the encounter with the other ones. The absence of any change in velocity, shape, or displacement 
evidences the no-scattering feature of these lumps, be they commutative or not. Summarizing, 
what survives at \t\ — > 00 is at most a one-soliton configuration with modified parameters and 
multiplied with constant unitary matrices. 

The asymptotic study just performed covers the £7(1) case as well, with some additional quali- 
fications. At rank one and m=2, the static lump is parametrized by 

\ Zu 0) = \ Zl ) = e 2iat ^ ia |0) , while 

|* 2 , i) = e Z2ct -* 2C U(t)\0) = C/(t)e 22at -^ a |0) (7.19) 

representing the moving lump remains time-dependent as \t\ — ► 00. Yet, the analysis is much 
simpler than for the nonabelian case. Since [39] 

{ Zl ,0\z 2 ,t) = e -^l 2 ' 2 (l-I|£ 2 | 2 ) , (7.20) 

the overlap s = (z\,0\z2,t) goes to zero exponentially fast and the solitons become well separated. 
Therefore, (7.7) tends to 

$t _> 1 _ 2\ Zl ){ Zl \ - p\z 2 ,t}(z 2 ,t\ , (7.21) 

which decomposes additively (rather than multiplicatively) into two independent solitons. 

Energies. It is technically difficult to compute the energy of true multi-solution configurations 
such as (7.4). However, because energy is a conserved quantity, we may evaluate it in the limit 
of \t\ — > 00. From the result (7.18) we seem to infer that generically only the £th soliton contributes 
to the total energy. Yet, this is incorrect, as can be seen already in the commutative case: Even 
though the moving lumps have disappeared from sight at \t\ = 00, their energies have not, because 
the integral of the energy density extends to spatial infinity, and so it is not legitimate to perform 
the large-time limit before the integration. What (7.18) shows, however, is the large-time vanishing 
of the overlap between the lumps in relative motion, so that asymptotically the lumps pertaining to 
different values of k are well separated. The total energy, being a local functional, then approaches 
a sum of m contributions, each of which stems from a group of r lumps, in isolation from the other 
groups. 

To compute the £th contribution, we must ignore the influence of the other groups of lumps. 
To this end we take the large-time limit and treat the £th group as a moving one-soliton. From 
(6.11) we know that its energy is 

Ee(v t ) = f(ve)E e (0) = 8irf(v e )q e . (7.22) 

By adding all contributions, it follows that the total energy of the m-soliton solution at any time is 

m 

E = 8tt qk cosh rjk sin ip^ ■ (7.23) 
fc=i 
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8 Conclusions 



In this paper we have demonstrated that the power of integrability can be extended to noncom- 
mutative field theories without problems. For the particular case of a 2+1 dimensional integrable 
sigma model which captures the BPS sector of the 2+1 dimensional Yang-Mills-Higgs system, the 
'dressing method' was applied to generate a wide class of multi-soliton solutions to the noncom- 
mutative field equations. We were able to transcend the limitations of standard noncommutative 
scalar and gauge theories, by providing a general construction scheme for and several examples 
of explicit analytical multi-soliton configurations, with full time dependence, for finite 6 values, 
and for any U(n) gauge group. Moreover, two proofs of large-time asymptotic factorization into 
a product of single solitons were presented, and the energies and topological charges have been 
computed. 

The model considered here does not stand alone, but is motivated by string theory. As was 
explained in two previous papers [41, 27], n coincident D2-branes in a 2+2 dimensional space-time 
give rise to open N=2 string dynamics [42 - 46] which (on tree level 2 ) is completely described 
by our integrable U(n) sigma model living on the brane. In this context, the multi-solitons are to 
be regarded as DO-branes moving inside the D2-brane(s). Switching on a constant NS 5-field and 
taking the Seiberg-Witten decoupling limit deforms the sigma model noncommutatively, admitting 
also regular abelian solitons. 

Since the massless mode of the open N=2 string in a space-time-filling brane parametrizes 
self-dual Yang- Mills in 2+2 dimensions [42], it is not surprising that the sigma-model field equa- 
tions also derive from the self-duality equations by dimensional reduction [28]. Moreover, it has 
been shown that most (if not all) integrable equations in three and less dimensions can be ob- 
tained from the self-dual Yang-Mills equations (or their hierarchy) by suitable reductions (see e.g. 
[49 - 53] and references therein). This implies, in particular, that open N=2 strings can provide 
a consistent quantization of integrable models in 1<D<3 dimensions. Adding a constant S-field 
background yields noncommutative self-dual Yang-Mills on the world volumes of coincident D3- 
branes and, hence, our noncommutative modified sigma model on the world volumes of coincident 
D2-branes [41]. It is certainly of interest to study also other reductions of the noncommutative 
self-dual Yang-Mills equations to noncommutative versions of KdV, nonlinear Schrodinger or other 
integrable equations, as have been considered e.g. in [54, 20, 22, 23]. 

Another feature of these models is a very rich symmetry structure. All symmetries of our 
equations can be derived (after noncommutative deformation) from symmetries of the self-dual 
Yang-Mills equations [55, 56] or from their stringy generalizations [38, 57]. For this one has to 
extend the relevant infinitesimal symmetries and integrable hierarchies to the noncommutative 
case. In fact, such symmetries can be realized as derivatives w.r.t. moduli parameters entering 
the solutions described in Sections 5 to 7 or, in a more general case, as vector fields on the moduli 
space. 

In this paper we considered the simultaneous and relative motion of several noncommutative 
lumps of energy. Like in the commutative limit, it turned out that no scattering occurs for the 
noncommutative solitons within the ansatz considered here. This means that the different lumps of 
energy move linearly in the spatial plane with transient shape deformations at mutual encounters 
but asymptotically emerge without time delay or changing profile or velocity. However, start- 
ing from the ansatz (4.12) containing second-order poles in £, one can actually construct explicit 

2 for a one- loop analysis see [47, 48] 
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solutions featuring nontrivial scattering of noncommutative solitons. Such configurations will gen- 
eralize earlier results obtained for commutative solitons [58 - 61] of the modified sigma model in 
2+1 dimensions. 

It would also be very interesting to construct the Seiberg-Witten map for our soliton configu- 
rations along the lines presented in [62]. This should allow one to compare their representation in 
commutative variables with their profile in the star-product formulation. 

Finally, an exciting relation between the algebraic structures of noncommutative soliton physics 
and string field theory has surfaced recently [63, 64, 65]. It is tempting to speculate that the 'sliver 
states' will play a role in the analysis of the open N=2 string field theory [66, 57]. We intend to 
investigate such questions in the near future. 
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A Coherent and squeezed states 

We collect some definitions and useful relations pertaining to coherent and squeezed states [39]. 
Coherent states for the Heisenberg algebra [a, a' ] = 1 are generated by the unitary shift operator 

D(J3) := e^-? a = e"^ ^ e~^ a . (A.l) 
It acts on the Heisenberg algebra as a shift, 




(A.2) 



and operates on the Fock space via 

£>(/?) |0) = e-5#V at |0) , (A.3) 

(0| J D t (/?i) J D(/3 2 )|0) = e-^e'^ 2 ^ 2 . (A.4) 

A representation of su(l, 1) can be constructed in terms of bilinears in the Heisenberg algebra, 

K + := \a) 2 , K_ := \a 2 , K := \{aa) + a) a) , (A.5) 

which fulfill 

[K ,K±] = ±K± and [K + ,K-] = -2K . (A.6) 
The unitary squeezing operator S may be defined as 

S(a) := e aK +-" K - = e^ K + e~ ln (i-ft)^o e -(K. = . s ^ ^ (A J) 

where 

a =: e y r and £ = e il? tanhr . (A.8) 
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S induces a representation on the Heisenberg algebra via 

„, 



, ,, . , ,, , / coshr — e^sinh-A 

S Tt K = f "t With 9 = . , , ■ (A.9) 



The action of S on the Fock space is characterized by 

5(0 io> = (i-ee) 1/4 e^+|o> , 

(0|5t(a)5(6)|0) = (i-e"i6) 1/4 (i-66) 1/4 (i-ei6r 1/2 • (A.10) 

The combined action of shifting and squeezing leads to squeezed states of the form 

I?,/?) := S(OD((3)\0) = e-l^a-eO^e^+^lO) , (A.ll) 

which are normalized to unity due to the unitarity of S and D. The corresponding transformation 
of the Heisenberg generators reads 



c \ „ ^ / a 



coshr — e 11? sinhr\ / a\ ( (3 



where the order of action was relevant. It is obvious that the unitary transformations defined 
in (4.26) are generated by Uk(t) = S(£k)D((3kt) with parameters {£,k,Pk} determined by 

B Asymptotic factorization for the two-soliton configuration 

We set out to proof the result (7.13). With (ii=—i and fj,2=^¥ : ~h we first consider the distinguished 
(z,z, t) coordinate frame. Let the entries of T 2 be polynomials of (the same) degree q 2 in c 2 - At 
large times, (4.26) implies that C2 — > — #2 and the i 92 term in each polynomial will dominate, 

T 2 -► t 92 f^l =: i 92 7 , (B.l) 

V 7n / 

where 7 is a fixed vector in group space. Due to the homogeneity of (7.4), the proportionality 
factor t q2 can be dropped and 7 be normalized to unity, 7^7 = 1. We may then substitute 

T\T 2 1 and P 2 -► T 2 7 7 + =: n = II* , (B.2) 

the projector onto the 7 direction. Note that the entries of the constant matrix II are commuting. 
Furthermore, we introduce p := 7 'Pi 7 which is a group scalar but noncommuting. Performing 
these substitutions in (7.4), the two-soliton configuration becomes 

*t _> 1 _ Mil Tl ^^ t} + ^oT x ^-^ rjn 

Mi T/(l-<rn)Ti M2 T/(l-<7n)Ti 

+ Bl a J_ U P 1 - ^^—U . (B.3) 
Hi 1—ap H2 l-crp 

In order to turn the terms in the first line to projectors, it is useful to introduce a constant matrix 

M = «n + (l-n) with a = ^—4 so that M f M = 1 - oTI (B.4) 

fj, — 1 
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and rewrite 

MTi =: f so that T x — f — T/ = — P^-r with P:=fi? . (B.5) 

t\m^MT x M Aft TtT 

We thus arrive at 

$t _> 1 _ ^-LpJL + ^i-pin + ^a-^—UPi - f^-^—U , (B.6) 
Hi M Mt n 2 M a m 1-ap fi 2 1-ap v ; 

after making use of 

/(M)n = n/(M) = /(a)n but /(M)(i-n) = (i-n)/(M) = l-n . (B.7) 

Further analysis is facilitated by decomposing <J?t [ n i Q parts orthogonal and parallel to 7. We insert 
our values (7.12) for the coefficients and find 

(i-n) $ f (i-n) -► (i-n) (1- 2 p) (i-n) (b.8) 



(i-n)*^ _ (i-n)(-2pi + ^^pi)7 

V a u(u+i) a) 



^t $ t (1 _ n) _> f tf_ 2 Ip + ^ 1 Pl ) ( i_n) 

V a fi— 1 1—ap J 



7^7 _ ^tf 1 _ 2 lpl + 2 H^)lpl + 2j^)^ 

V a a /i(/i+i) a a /i— i 1— op p 1—ap/ 

With the help of the identities 

7 t P 7 = a 7 t — — Pi j a = a a — - — and 

1 — crPin 1 — ap 

7tp(i_n) = aft- ^— -Pi(l-n) = a — ^7 f A(l-n) (B.9) 

1 — a/ill 1 — 07? 

plus some algebra, the pieces of <J>t m (B.8) simplify to 

(i-n) *t (i-n) -► (i-n) (1-2P) (i-n) 
(i-n)$t^ _^ (i-n) (-2P + 2pP)7 
7t$t(!_n) ft (-2P) (i-n) 

^t$t^ _^ (1 _ p _ 2P + 2pP) 7 , (B.10) 

which proves the factorization 3 

$t (i-2P)(i-pn) . (B.11) 

At |t| — ► 00 we are thus left with a static one-soliton configuration with modified parameters 
(Ti — > MTi) and multiplied with a constant matrix (1— /?n). If we move with the second lump the 
proof is the same, because then a — ► and Pi (in place of P2) goes to a constant projector. 

3 We have omitted writing 'soliton indices' for P and IT, to avoid cluttering our formulae. 
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C Asymptotic form of m-soliton configurations 

The dressing method may be used to show the large-time factorization of multi-soliton solutions 
into a product of one-soliton solutions directly on the level of the auxiliary function ip(t,x,y,Q. 

Starting from the one-soliton ansatz, 

* = 1 + £=^A (CI) 

with a hermitian projector Pi, we may try to obtain a two-soliton solution by 'dressing' the latter 
with another factor of the same kind / fj, 2 ), 

^ = + + • (C2) 

The reality condition (4.2) is identically fulfilled as soon as P\ and P2 are hermitian projectors, 

Pi = Ti^U^i and h = T 2 — l^lj . (C.3) 

T/Ti T 2 f r 2 

The removability of the pole for £ — ► Jx 2 in (4.10) and (4.11) is guaranteed if 

(1-^2)02^2 = =► c 2 f 2 = f 2 ^ 2 , (C.4) 

qualifying the second factor on the r.h.s. of (C.2) as a standard one-soliton solution, as obtained 
from solving (4.31). We therefore redenote P2 — > P 2 and T 2 — > T 2 . 4 The vanishing of the £ — ► /ii 
residue, on the other hand, leads to 

(l-A)(l + ^ftW(l + ^£)A = , (C.5) 

which does not yield a 'holomorphicity condition' for T\ but rather demonstrates that P\ is not a 
standard one-soliton projector. 

Let us move with the first lump. In the \t\ — > 00 limit, then, the arguments used in (B.l) 
and (B.2) apply, meaning that P 2 — > n 2 . Since II 2 is a coordinate-independent projector, the ci 
in (C.5) can be moved next to P\ which implies 

p 2 ^n 2 (l-Pi)ciPi = = f 1 z 1 , (c.6) 

so that the large-time limit of Pi is also a standard one-soliton projector, and we again redenote 
Pi — > Pi. Hence, 

iim^ 2 = (i + aza^)(i + ^n 2 ) , (C.7) 

|t|-oo V C-W ' V / 

which yields 

lim $ 2 = (l-piPi) (l-p 2 n 2 ) . (C.8) 

|t|— >oo V / \ / 

On the other hand, if we move with the second lump, then, in the large-time limit, c\ — > —((3i—(3 2 )t, 
implying Pi — > LTi where 111 is a coordinate-independent projector. In this limit, (C.5) is satisfied, 
and (C.4) shows that P 2 =P 2 describes a one-soliton configuration. Hence, in this frame we find 

lim $t = (l-^ni) (l-p 2 P 2 ) . (C.9) 

\t\— >oo V / \ / 



It will be seen shortly that these objects coincide with the ones defined in Appendix B. 
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It can be shown, in fact, that the multi-soliton ansatz (4.13) is equivalent to the product ansatz 

m _ ^ 

= ~V\\l + ^^P k ] , with P k = f k ^f\ , (CIO) 
11 L C-Uh J T]T h 



k=i s w J 1 k ± k 

where in general the are 'non-holomorphic' matrices. By induction of the above argument one 
easily arrives at the m-soliton generalization of (C.8). Namely, in the frame moving with the £th 
lump we have 

lim & m = (1 - Pl U 1 ) . . . (1 - pt-iUt-!) (1 - piPt) (1 - Pi+lUt+l) ■ ■ ■ (1 - Pm^m) ■ (C.ll) 

\t\— »oo 

This provides an alternative proof of large-time factorization. 
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